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Abstract: 

Steady  plastic  shocks  generated  by  planar  impact  on  metal-polymer  laminate  composites,  are 
analyzed  in  the  framework  of  gradient  plasticity  theories.  The  laminate  material  has  a  periodic 
structure  with  unit  cell  composed  of  two  layers  of  different  materials.  First  and  second  order 
gradient  plasticity  theories  are  used  to  model  the  structure  of  steady  plastic  shocks.  In  both 
theories,  the  effect  of  the  internal  structure  is  accounted  for  at  the  macroscopic  level  by  two 
material  parameters  depending  upon  the  layer’s  thickness  and  the  properties  of  constituents. 
Those  two  structure-parameters  are  shown  to  be  uniquely  determined  from  experimental  data. 
Theoretical  predictions  are  compared  with  experiments  for  different  cell  sizes  and  for  various 
shock  intensities.  In  particular,  the  following  experimental  features  are  well  reproduced  by  the 
modeling: 

the  shock  width  is  proportional  to  the  cell  size; 

the  magnitude  of  strain  rate  is  inversely  proportional  to  cell  size  and  increases  with  the 
amplitude  of  applied  stress  following  a  power  law. 

While  these  results  are  equally  described  by  both  the  plasticity  theories,  the  first  gradient 
plasticity  approach  seems  to  be  favored  when  comparing  the  structure  of  the  shock  front  to 
experimental  data. 


1)  Introduction 

The  visco-plastic  response  of  metals  and  polymers  at  high  loading  rates  can  be  investigated 
using  shock  wave  experiments.  Plate  impact  experiments  have  provided  a  large  quantity  of 
experimental  data  on  steady  plastic  shocks,  especially  for  polycrystalline  metals,  Rice  et  al 
(1958),  Johnson  and  Barker  (1969),  Swegle  and  Grady  (1985).  By  analyzing  a  stationary  plastic 
shock,  the  following  data  can  be  obtained:  (i)  stress  amplitude,  (ii)  particle  velocity  versus  time 

at  a  given  point,  (iii)  plastic-shock  velocity  C,  (iv)  particle  velocity  behind  the  shock  front 

and  (v)  maximum  strain  rate  within  the  shock  front  layer.  The  Hugoniot  curve  is  an  important 
material  characteristic  obtained  from  shock  analysis.  From  data  (iii)  and  (iv),  a  linear  relationship 

between  C  and  v  is  obtained  from  which  the  Hugoniot  can  be  deduced  (see  for  example, 

Molinari  and  Ravichandran  (2004)). 

The  problem  of  relating  shock  experimental  data  to  material  properties  is  a  problem  of 
fundamental  importance  which  was  studied  in  detail  for  certain  classes  of  materials;  see  for 
example,  Clifton  (1971),  Swegle  and  Grady  (1985),  Molinari  and  Ravichandran  (2004).  The 
shock  front  is  a  thin  layer  where  particle  velocity,  stress  and  strain  sustain  a  rapid  variation.  The 
spatial  distribution  of  those  quantities  within  the  shock  front,  which  is  usually  referred  to  as  the 
shock  structure,  is  in  many  cases  controlled  by  viscous  effects.  For  polycrystalline  metals, 
viscous  effects  can  be  related  to  dislocation  gliding.  For  more  details  and  references,  it  is  referred 
to  Molinari  and  Ravichandran  (2004),  where  a  detailed  analysis  is  given  for  moderate  shocks  in 
metals,  with  a  specific  application  to  aluminum.  As  mentioned  before,  the  maximum  of  the  strain 
rate  within  the  plastic  shock  layer  can  be  measured  for  a  given  stress  amplitude  A  a  (stress 

jump  across  the  shock  front).  A  particular  definition  of  the  maximum  strain  rate  ( sSG )  was 

proposed  by  Swegle  and  Grady  (1985)  as  being  the  ratio  of  the  maximum  time  derivative  of  the 
particle  velocity  to  the  steady  shock  velocity.  Experimental  data  were  found  by  Swegle  and 

Grady  (1985)  to  follow  closely  the  relationship  sSG  =  B(Aa)hsG  with  the  exponent  hSG  =  4 
being  nearly  the  same  for  all  materials  considered.  These  authors  related  hSG  to  the  material 

strain  rate  sensitivity  m  through  the  relationship,  hSG  =  2/ m .  Thus,  it  was  concluded  that 

materials  exhibit  nearly  the  same  strain  rate  sensitivity  m  =  0.5  for  the  deformation  (strain  rate) 
regime  encountered  in  the  shock  experiments  considered.  Molinari  and  Ravichandran  (2004) 
have  obtained  an  analytical  solution  for  steady  state  shocks  of  moderate  intensity  in  metals.  From 


their  analysis  it  was  found  that  the  maximum  strain  rate  within  the  shock  can  be  related  to  the 

stress  jump  by  a  power  law  |£|  =  B{ Act)*  with  h  =  1 .87(1  Jm)  +  0.62  for  aluminum.  Therefore, 

for  aluminum,  the  value  h  =  4  providing  a  good  fit  of  experimental  data,  according  to  the 
analysis  of  Molinari  and  Ravichandran  (2004),  to  the  value  m=  0.561  of  the  strain  rate 
sensitivity  for  the  high  loading  rates  generated  in  shock  experiments.  Therefore,  having  a  unique 
value  of  the  strain  rate  sensitivity  (m  =  0.5  according  to  Swegle  and  Grady  (1985))  for  all 
metallic  materials,  is  debatable. 

In  any  case,  those  findings  are  related  to  the  particular  form  of  the  constitutive  law  chosen  for 
the  modeling  and  more  particularly  to  the  specific  form  of  the  “nonlinear-viscous”  response 
within  the  shock  layer.  Indeed,  it  is  frequently  assumed  that  the  shock  structure  is  controlled  by 
the  dissipative  effects  associated  to  the  material  constitutive  response.  However,  it  is  well  known 
that  shocks  can  also  be  structured  by  dispersive  effects.  A  good  illustration  is  given  by  the 
Korteweg-Devries  equation  governing  shallow  water  waves.  In  general,  dispersive  effects  result 
from  wave  interaction  with, 

(i)  external  boundaries  of  the  body  (e.g.  dispersion  of  elastic  waves  in  a  circular  rod, 

Pochhammer  (1876)  and  Chree  (1889)  ); 

(ii)  inner  boundaries  due  to  material  heterogeneities  (inclusions,  layers). 

In  this  paper,  the  primary  faocus  is  on  plane  waves  propagating  in  unbounded  media,  and 
therefore  interactions  of  type  (i)  shall  not  be  considered. 

In  the  last  four  decades,  nonlocal  linear  elastic  laws  (in  particular  gradient  deformation 
theories)  were  formulated  with  the  purpose  of  describing,  within  a  phenomenological  approach, 
wave  dispersion  related  to  the  internal  structure  of  heterogeneous  materials,  Green  and  Rivlin 
(1964),  Eringen  (1972,  1976),  Mindlin  (1964),  Kunin  (1982,1983),  Luciano  and  Willis  (2001). 

For  particular  composite  materials  such  as  elastic  multilayer  periodic  media  or  fiber 
composites,  wave  dispersion  can  be  directly  analyzed  without  using  macroscopic  nonlocal 
formulations.  By  considering  the  basic  principles  of  continuum  mechanics,  wave  dispersion  of 
hannonic  waves  can  be  treated  by  analytical  means.  From  the  dispersion  curve,  the  effects  of 
elastic  properties  of  constituents  can  be  characterized  as  well  as  those  of  the  geometry  of  the 
internal  structure  (e.g.,  layer’s  thickness  of  laminates),  Sun  et  al.  (1968),  Achenbach  (1973), 
Christensen  (1979),  Boutin  (1996). 

Analyses  of  wave  dispersion  in  nonlinear  materials  accounting  for  the  effects  of  internal 
structure  are  less  advanced.  Recently,  Grady  (1998)  analyzed  steady  shock  waves  in 
polycrystalline  materials  and  considered  the  possibility  for  shocks  to  be  structured  by  dispersive 
effects  due  to  wave  scattering  on  grain  boundaries.  These  dispersive  effects  were  viewed  as  an 
alternative  to  visco-plasticity  for  structuring  steady  shock  profiles.  A  comprehensive  review  of 
experimental  work  and  numerical  analysis  of  shocks  in  laminated  materials  can  be  found  in 
Nesterenko  (2001),  Chapter  3,  with  a  detailed  discussion  of  transient  shock  waves.  It  is  shown 
how  the  leading  wave  is  attenuated  by  wave  reflection  at  the  layer’s  interfaces.  However,  as  for 
the  Korteweg-Devries  equation  is  considered,  the  main  signal  appears  finally  to  be  the  result  of 
the  competition  between  non-linear  effects  and  dispersion  effects,  Whitham  (1974). 


The  aim  of  the  present  work  is  to  analyze  the  effects  of  wave  dispersion  on  steady  shock 
propagating  in  layered  periodic  media  with  nonlinear  elastic  or  plastic  properties.  The  direction 
of  propagation  is  assumed  to  be  perpendicular  to  the  layers.  Dispersive  effects  are  due  to  wave 
reflection  and  refraction  at  the  layer's  interfaces.  A  detailed  analysis  of  these  interactions,  which 
was  feasible  for  linear  elastic  materials,  is  hard  to  achieve  for  nonlinear  responses  by  pure 
analytical  means,  but  insight  can  be  gained  through  numerical  simulations,  Benson  and 
Nesterenko  (2001).  The  framework  adopted  here  is  based  on  a  simplified  constitutive  assumption 
allowing  exploring  by  simple  calculations  the  structure  of  plastic  steady  shocks.  It  is  assumed 
that  the  effect  of  the  layered  structure  can  be  accounted  for  at  the  macroscopic  level  by  a 
non-local  constitutive  model. 

Many  efforts  have  been  recently  devoted  to  the  development  of  gradient  plasticity  theories. 
These  approaches  are  aimed  at  accounting  for  the  effects  of  the  internal  structure  on  the  overall 
effective  mechanical  response,  Aifantis  (2003),  Fleck  and  Hutchinson  (1993,  2001),  Gao  et  al. 
(1999),  Huang  et  al.  (2000).  Among  the  different  formulations  developed  in  the  literature,  two 
gradient  plasticity  approaches  are  considered  and  their  applicability  in  modeling  the  propagation 
of  steady  shock  waves  in  laminates  is  evaluated. 

The  modeling  of  shock  wave  experiments,  recently  performed  by  Zhuang  et  al.  (2002,  2003) 
on  laminated  materials,  is  considered  here.  In  these  experiments,  the  stress  jump  A  a  across  the 

shock  and  the  maximum  strain  rate  within  the  shock  |Z|  ,  are  found  to  be  related  by  a  power 

law,  If]  =  B(Aaf' ,  as  for  polycrystalline  metals.  However,  the  value  of  the  exponent  is  now 

in  the  range  1.8  <hL  <2.4  depending  on  the  specific  laminate  considered.  These  results  are 

noticeably  different  from  the  value  h  =  4  found  in  shock  experiments  on  single  phase 
polycrystalline  metals  by  Swegle  and  Grady  (1985).  The  modeling  proposed  in  this  paper 

provides  a  quite  good  representation  of  the  experimental  results  of  Zhuang  et  al  (2002,  2003). 

2)  Plastic  shocks  in  laminates 

Shock  waves  generated  in  laminate  composites  by  planar  impact,  have  been  studied 
experimentally,  see  for  example,  Lundergan  and  Drumheller  (1971)  and  Zhuang  et  al  (2002, 
2003).  These  waves  propagate  in  the  impact  direction  perpendicular  to  layers  as  shown  in  Fig.  1. 
The  laminates  considered  in  this  paper  are  constituted  by  alternate  layers  of  materials  (1)  and  (2) 

of  thickness  Li  and  Z?  respectively.  The  elementary  cell,  of  total  thickness,  L=Ll+L2 ,  is  made  of 

the  two  layers  of  dissimilar  materials. 

Beyond  a  certain  propagation  distance,  quasi-steady  plastic  shocks  are  fonned  (the  meaning 
of  “quasi-steady”  will  be  discussed  later).  The  shock  structure  should  depend  on  the  properties  of 
the  constituents  and  also  on  the  layer's  thickness.  Dispersion  effects  due  to  wave  reflection  and 
refraction  at  the  layer’s  interfaces,  are  thought  to  play  a  significant  role  in  structuring  the  shock 


front  together  with  nonlinear  effects. 

The  propagation  distance  for  having  steady  plastic  shocks  depends  intimately  on  the 
properties  of  the  laminate  constituents.  In  the  experiments  by  Lundergan  and  Drumheller  (1971) 
on  steel-PMMA  laminates,  it  was  observed  that,  after  a  propagation  distance  2L  of  two  cell  sizes, 
the  the  particle  velocity  versus  time  was  almost  identical  at.  In  general,  at  two  homologous 
points  distant  by  a  multiple  of  L  (periodicity  of  the  internal  structure)  recordings  of  the  particle 
velocity  should  be  identical  (up  to  a  time  shift).  Observations  similar  to  those  of  Lundergan  and 
Drumheller  (1971)  were  reported  by  Zhuang  et  al  (2002,  2003)  in  their  study  of 
polycarbonate-steel  laminates.  In  the  experiments  made  by  Nesterenko  et  al.  (1983)  on 
aluminum-copper  laminates,  it  was  shown  that  after  a  distance  of  6L  no  stationary  wave  profile 
was  found.  The  higher  impendence  mismatch  between  constituents  realized  in  metal-polymer 
laminates  is  probably  a  reason  for  having  stationary  shocks  after  a  short  traveling  distance. 

Zhuang  et  al.  (2002,  2003)  have  measured  the  stationary  propagation  speed  C  of  the  main 
shock  in  polycarbonate-steel  laminates.  Let  us  consider  now  an  observer  moving  at  the  constant 
velocity  C.  For  this  observer,  the  wave  profile  (spatial  distribution  of  the  particle  velocity) 
fluctuates  with  a  time  periodicity  equal  to  L/C,  and  a  steady  wave  profile  can  be  defined  by  time 
averaging  over  the  period  L/C  of  the  time-fluctuating  profile. 

The  aim  of  this  paper  is  to  analyze  these  steady  shock  fronts.  The  details  of  wave  reflection 
and  refraction  at  layer's  interfaces  are  not  analyzed.  Rather  the  steady  shock  defined  above  by 
time  averaging,  is  supposed  to  be  representative  of  the  dynamic  response  of  the  laminate 
structure  after  a  certain  propagation  distance.  This  steady  shock  moves  with  the  constant 
velocity  C  and  maintains  a  constant  shape. 

Assuming  periodicity,  a  steady  profile  can  be  defined  in  a  different  way.  For  the  observer 
moving  with  the  constant  speed  C,  material  properties  are  seen  in  average.  It  is  thought  that  the 
average  material  response  experienced  by  the  observer  can  be  identified  with  the  macroscopic 
effective  response  of  the  laminate.  In  this  paper,  a  phenomenological  gradient  plasticity 
approach  is  used  to  describe  the  effective  response  of  the  laminate,  and  steady  shock  waves  are 
analyzed  in  this  framework. 

The  constant  wave  profile  in  a  steady  shock  is  the  result  of  the  balance  between  steepening 
effects  related  to  the  non-linearity  of  the  material  response  and  smoothing  of  the  wave  front  due 
to  dissipative  and/or  dispersion  effects.  The  dispersion  effects  are  expected  to  be  important  in 
case  of  a  large  contrast  between  the  mechanical  properties  of  constituents. 

For  large  impact  velocities,  the  material  undergoes  compressive  stresses  beyond  the  initial 
yield  limit  of  the  constituents  and  a  plastic  shock  wave  is  formed.  Plastic  shock  waves  have  been 
studied  in  the  experimental  work  of  Zhuang  et  al.  (2002,  2003).  For  a  given  shock  strength 
(stress  amplitude),  the  particle  velocity  was  recorded  in  terms  of  time  at  different  locations  and 
the  velocity  C  of  the  stationary  plastic  shock  was  measured.  It  was  shown  that  the  shock  velocity 

C  varies  linearly  with  the  particle  velocity  v  measured  at  the  rear  of  the  plastic  shock  front.  In 
addition,  the  maximum  strain  rate  (maxUI)  within  the  shock  layer  was  measured  and  the 


following  power  law  was  obtained, 

max|f|  =  b(ct  (1) 

where  cr  is  the  compressive  stress  at  the  rear  of  the  plastic  shock.  The  exponent  was  shown  to 

be  in  the  range,  1 .5  <  hL  <  2.2 ,  for  the  different  laminates  tested. 

The  power  law  (1)  for  a  two  phase  laminate  is  similar  to  that  of  a  single  phase  metal  which, 
according  to  Swegle  and  Grady  (1985),  can  be  written  as, 

max^l  =  B(Acr)''  (2) 

where,  Act  =  <j ~  -  a+  is  the  stress  jump  across  the  plastic  shock  ( <j+  is  the  compression  stress 
ahead  of  the  plastic  shock).  Note  that  for  the  polymer-metal  laminates  considered  here,  <j+  is 
small  or  zero  and  thus  A  o'  *  a  .  For  single  phase  metals,  the  value  of  the  exponent  h  in  the 

power  law  (2)  was  found  to  be  close  to  h  «  4  which  is  noticeably  larger  than  values  of  hL  for 

laminates.  Dispersion  effects  due  to  wave  reflection  and  refraction  at  the  laminate  material 
interfaces  may  explain  this  difference.  Indeed,  the  scope  of  this  paper  is  to  analyze  how  the 
shock  profile  can  be  affected  by  dispersion  effects.  In  an  initial  analysis,  dissipative  effects  are 
neglected  so  as  to  analyze  the  role  of  dispersion  separately. 

Laminar  composites  are  considered,  which  are  made  of  two  dissimilar  materials  with  “soft” 
(e.g.,  polymer  such  as  polycarbonate  (PC))  and  “hard”  (e.g.,  metal  such  as  stainless  steel  (SS)) 
layers  with  volume  fractions  being  respectively, 

fi=L2IL .  (3) 

Material  characteristics  of  the  constituents  considered  in  the  modeling  are  given  in  Table  1 . 

Two  different  composites  are  obtained  by  varying  the  layer’s  thickness.  The  PC37/SS19 

laminate  corresponds  to  Lx  =  0.37  mm  and  L2  =0.1 9  mm  and  PC74/SS37  is  obtained  by 

doubling  the  width  of  layers. 

Experimental  data  (Zhuang  et  al.  (2002,  2003))  plotted  in  Fig.  2  show  that  the  Lagrangian 
shock  speed  is  related  to  the  particle  velocity  v  at  the  rear  of  the  shock  by  a  linear  relationship. 
For  plastic  shocks  of  moderate  amplitude,  this  linear  relationship  can  be  written  as: 

C  =  c+  +  S(v~  -v+) 


(4) 


where  the  subscript  (+)  refers  to  the  state  ahead  of  the  plastic  shock;  v+is  the  particle  velocity 

ahead  of  the  shock  front  and  c+  is  the  sound  velocity  corresponding  to  state  (+). 

Note  that  the  experimental  data  shown  in  Fig.  2  provide  identical  linear  relationships  for 
both  composites,  showing  no  sensitivity  to  the  cell  size.  The  sensitivity  of  the  shock  wave  profile 
to  the  cell  size  is  analyzed  in  Nesterenko  (2001).  For  solid  teflon-paraffin  laminates,  Nesterenko 
et  al.  (1983)  report  that  the  amplitude  of  the  particle  velocity,  at  a  distance  of  2L  or  4L  from  the 
impact  plane,  is  nearly  insensitive  to  the  cell  size.  It  is  believed  that  this  independence  with 
respect  to  cell  size,  as  in  the  experiments  of  Zhuang  et  al.  (2002,  2003),  is  a  feature  of  steady 
shocks.  For  transient  shock  profiles  there  is  a  strong  effect  of  the  cell  size  on  the  amplitude  of  the 
leading  signal,  see  Nesterenko  (2001)  Fig.  3.6  for  aluminum-copper  laminate. 

The  slope  of  the  best  fit  lines  for  the  experimental  data  in  Fig.  2  is  given  by, 

S  =  2.l  (5) 

The  relationship  (4)  with  S=2.\,v+  =  18.2  mis  (value  determined  later)  and  c'  =1.600 /h/.v 

shows  a  good  fit  with  the  experimental  data.  Note  that  the  Hugoniot  curve  of  the  laminate  can  be 
obtained  from  (4),  for  further  discussion,  see  section  4. 


3)  Basic  equations 

A  Lagrangian  formulation  is  used  to  account  for  the  large  uniaxial  defonnations  occurring  in 
shock  wave  experiments.  The  direction  of  propagation  of  the  plane  wave  is  denoted  by  Ox,  Fig.  1; 
u  is  the  component  of  the  displacement  in  this  direction,  other  components  along  directions  y  and 

z  being  equal  to  zero.  The  orthogonal  frame  Oxyz  has  unit  vectors  e, ,  e2  and  e3 .  Note  that  the 

stresses  and  defonnations  considered  in  the  following  are  defined  at  the  macroscopic  level  of  the 
composite. 

The  present  position  x  of  a  particle  is  given  by: 

x  =  X  +  u(X)  (6) 


where  X  is  the  initial  position  prior  to  any  defonnation.  The  deformation  gradient  has  the  form, 


F_  =  Alel  0  ex  +  e2  0  e2  +  e3  0  e3 


where  the  longitudinal  stretch  is, 

,  dx  ,  du 
A,  =  —  =  1  +  — 


dX 


dX 


(V) 

(8) 


The  standard  multiplicative  decomposition  of  F_  into  elastic  (e)  and  plastic  parts  (p) ,  see  Fig.  3, 
can  be  written  as, 

F  =  Fe  Fp 


(9) 


with 


F_e =  ^ifi  ®  §_\  +  Ae2e2  ®  e2  +^2^3  ®^3  (10) 

Fp  =  Apex®ex+Ape2®e2+Ape3®e3  (11) 

Note  that  X2  =  A3  and  Ap  =  Ap  due  to  the  rotational  symmetry  around  the  direction  ex . 

The  condition  of  plastic  incompressibility: 

Ap(Ap)2=l  (12) 

and  the  multiplicative  decomposition  (9)  provide  three  independent  relations  relating  the 
stretches  Ax,  Xx,  X2 ,  Ap  and  Ap .  Therefore  all  stretches  can  be  expressed  in  terms  of  2, 

and/lf , 

Xx=AjAp,  X2  =  (2f)1/2 ,  Ap=(Apy1/2.  (13) 

Denoting  by  T  the  nominal  stress  tensor  with  respect  to  the  initial  configuration  C0 ,  one 
has,  Clifton  (1971),  Lubarda  (2002): 

T  =  (F_py\Te  (14) 


The  elastic  constitutive  law  is  expressed  in  the  intermediate  configuration  Cr ,  Fig.  3,  as: 


rj^e 


dyX_ 

dFe 


(15) 


where  t//'  is  the  macroscopic  elastic  strain  energy  (volume  average  of  the  layer’s  strain  energy). 


From  (11)  and  (14), 


(16) 


To  represent  the  nonlinear  elastic  response  for  compressive  loading,  the  elastic  law  (15)  should 
include  higher  order  elastic  coefficients.  Second  order  elastic  coefficients  can  be  characterized  by 
ultrasonic  measurements;  however  this  information  is  not  available  for  the  laminates  studied  here. 
Therefore  second  order  elastic  constants  will  be  introduced  as  parameters,  and  a  parametric 
analysis  will  be  perfonned  to  quantify  how  results  are  affected  by  second  order  elastic  constants. 

The  following  nonlinear  constitutive  law  is  adopted, 


Tx\  =  (Cxxxx£ex  +  2CXU2£2)(1  -  b{£[  +  0.2ffa')) 


(17) 


where  Cnn  and  Cn22  are  components  of  the  overall  (effective)  linear  elastic  moduli  of  the 

laminate.  For  both  laminates  PC37/SS19  and  PC74/SS37,  one  has  the  same  values: 

Cnu  =  6.84  GPa  ,  CU22  =  3.60  GPa  .  (18) 

Second  order  elastic  properties  are  scaled  by,  6=66.7  .  This  value  will  be  varied  in  the 
parametric  analysis. 

Elastic  deformations  e[  and  s\  are  expressed  in  terms  of  \  and  /If  by, 

<=log(^)  =  log(V^) 

si= log® = began  (is) 

Combining  (16),  (17)  and  (19),  the  stress  component  Tn,  denoted  henceforth  for  simplicity  as 
T ,  can  be  written  in  terms  of  and  /If , 

T  =  T(\,%). 

The  compressive  stress  (with  positive  sign)  is  defined  as  follows, 
a  =  -T . 

4)  Steady  plastic  shocks  and  the  Hugoniot 

Plate  impact  experiments  are  consdiered,  where  stationary  plastic  shocks  are  fonned  beyond 
a  certain  propagation  distance.  These  shocks  move  with  constant  speed  and  constant  shape.  Two 
cases  are  considered, 

(i)  the  stress  amplitude  exceeds  the  elastic  limit  of  the  composite  so  as  to  fonn  a  plastic  shock 
preceded  by  an  elastic  precursor  with  higher  propagation  speed; 

(ii)  the  stress  amplitude  is  sufficiently  large  to  form  a  single  elasto-plastic  shock. 

The  analysis  is  restricted  to  weak  shocks,  i.e.,  temperature  effects  are  neglected  and  the 
process  is  considered  as  isentropic,  see  Clifton  (1971)  and  Molinari  and  Ravichandran  (2004)  for 
further  comments. 

4.1)  Plastic  shock  with  elastic  precursor 

For  moderate  impact  velocities,  the  plastic  shock  is  preceded  by  an  elastic  precursor  where 
the  material  is  compressed  from  the  initial  state  (0)  to  state  (+)  corresponding  to  first  plastic 
yielding.  The  values  of  stress,  particle  velocity  and  stretches  associated  to  state  (+)  are  denoted 

by cr+ = -T+,v+,  X  and  Ap'  ,  where  the  subscript  (1)  has  been  dropped  to  simplify  the 
notation.  In  the  plastic  shock,  the  material  is  further  compressed  to  state  (-)  with  stress 


(20) 

(21) 


a  =—T  reached  at  the  rear  of  the  plastic  front,  see  Fig.  4.  The  transition  between  states  (+) 


and  (-)  occurs  within  the  plastic  shock  front,  a  thin  layer  where  the  material  sustains  high  strain 
rates. 

The  equation  of  conservation  of  momentum  is, 


dT 


dv 


- =  Pa  - 

dX  0  dt 


(22) 


where  t  is  time,  v  =  du(X,t)  /  dt  is  the  particle  velocity  and  pQ  is  the  mass  density  in  the 


reference  configuration,  C0 .  The  compatibility  equation  has  the  form, 

dv  _  dA 
dX  dt 

For  an  observer  moving  with  the  plastic  shock  velocity  C,  the  problem  is  steady  and  all 
variables  can  be  expressed  in  terms  of  the  parameter, 


(23) 


4  =  X-Ct  (24) 

with  states  (+)  and  (-)  corresponding  respectively  to  £  =  +oo  and  g  =  -oo . 

Using  (24),  the  equations  (22)  and  (23)  are  written  as, 
dT  dv 

Tr~Cp^  (25) 

^  =  (26) 

d%  dt, 

Upon  integration,  it  can  be  shown  that  any  state  within  the  plastic  shock  layer  has  to  satisfy  the 
following  relationships, 


T-T+=-p0C(v-v+ ) 


(27) 


v-v+  =-C(A-A+)  (28) 

By  combining  (27)  and  (28)  the  equation  defining  the  Rayleigh  line  is  obtained  (see  Fig.  4), 
T-T+  =p0C\A-A+)  (29) 

In  particular,  state  (-)  is  related  to  state  (+)  by, 

T--T+=-p0C(v--v+ )  (30) 

-  +=-C(/T-/ 1+) 


v  —v 


(31) 


T  -T+  =  p0C2(A-  -  A+) . 


(32) 


The  Hugoniot  is  the  curve  defined  by  the  ensemble  of  points  (a  ,A  )  obtained  by  varying 
the  stress  amplitude  cr  ,  Fig.  4.  The  intercept  of  the  elastic-loading  curve  A  —>  T(A  ,AP  =  1) 

with  the  Hugoniot  is  the  state  (+),  see  Fig.  4.  Thus,  the  compressive  stress  <r+  at  state  (+)  follows 
the  following  relationship, 

a+=-T(A+,l).  (33) 

The  state  (+)  is  also  defined  as  being  the  initial  yield  limit  of  the  composite.  It  will  be  shown  that 
the  precise  characterization  of  state  (+)  is  not  a  crucial  point.  This  is  related  to  the  fact  that,  for 
the  metal-polymer  composite  studied  here,  the  value  of  the  initial  yield  limit  is  small  and  does 
not  affect  the  results.  Therefore,  a  conventional  definition  of  state  (+)  can  be  adopted  by 

setting  A+  =  0.99 ,  see  Fig.  4.  The  corresponding  stress  is  cr+=0. 11 5GPa  .  The  particle  velocity  at 

state  (+)  is  given  by  (see  for  instance  Molinari  and  Ravichandran  (2004)), 

v+=cj  1-A+)  (34) 

where 


is  the  velocity  of  the  elastic  precursor.  p0  is  the  average  mass  density  of  the  composite  in  the 
reference  configuration: 

A  =/  P\  +fi  Pi  ■  (36) 

The  two  PC-SS  laminates  differ  only  by  a  factor  two  in  the  layer's  thickness;  therefore  they  have 
same  mass  density  /i0  =3,468  kg/m3 .  The  elastic  compression  curve  is  also  identical  and  the 


initial  yield  limit  is  reached  for  the  same  value  of  the  nominal  stress  denoted  byT+.  Both 


composites  have  the  same  state  (+).  The  particle  velocity  at  the  rear  of  the  elastic  shock  is 
v+  =18.2  ml  s  .  The  velocity  of  the  elastic  shock  precursor  isce/  =  1,820  ml  s  . 

By  substituting  (31)  into  (4), 


C  = 


1  -S(A+  -A-) 


(37) 


<7  =<7H(A)-  (4°) 

The  configuration  of  a  plastic  shock  with  an  elastic  precursor  is  physically  realizable  if  the 
velocity  ce/  of  the  elastic  precursor  is  larger  than  the  plastic  shock  velocity  C, 


The  case  where  this  condition  is  not  satisfied  (shocks  of  larger  amplitude)  is  studied  in  the  next 
section. 

Note  that  the  analytical  form  (39)  of  the  Hugoniot  was  derived  by  considering  shocks  of 
weak  enough  stress  amplitude.  It  is  assumed  in  the  following  that  the  Hugoniot  can  be  described 
in  first  approximation  by  the  same  expression  (39)  extrapolated  to  a  larger  range  of  stress 
amplitude.  The  consistency  of  this  hypothesis  will  be  validated  by  comparing  the  prediction  of 
the  theory  with  a  different  set  of  experimental  data. 


4.2)  Case  of  a  single  elastic-plastic  shock 


When  the  condition  (41)  is  not  satisfied,  an  unique  shock  is  formed,  relating  the  initial  state 
(0)  to  the  state  (-)  as  shown  in  Fig.  6  for  the  stress  amplitude,  cr  =  2 GPa  .  The  Rayleigh  line  is 
given  by  (29)  where  the  state  (+)  is  replaced  by  state  (0), 

T  =  p0C\A- 1).  (42) 

The  intercept  of  the  Rayleigh  line  with  the  elastic  loading  path  A  — »  a  =  -T(A,AP  =  1)  defines 

the  state  (*).  States  (0)  and  (*)  are  related  by  a  pure  elastic  path.  Plastic  flow  occurs  between 
states  (*)  and  (-). 

The  shock  velocity  is  given  by  (42), 


Figure  2  shows  the  shock  velocity  C  in  terms  of  the  particle  velocity  v  at  the  rear  of  the 

shock.  Experimental  data  are  represented  by  stars  for  PC37/SS19  and  by  squares  for  PC74/SS37. 
The  modeling  prediction  corresponds  to  the  solid  line  with  dots  representing  different  levels  of 

the  applied  stress  cr  (  from  1  GPa  to  9  GPa).  For  given  values  ofcr  ,  the  theoretical  results  have 
been  obtained  by  calculating  X  with  (39)  and  (40).  The  values  of  material  characteristics  in  (39) 
are  given  in  Table  2.  The  shock  velocity  C  is  obtained  with  (43)  and  with  (3 1)  where  state  (+) 
is  replaced  by  the  initial  state  (0),  v  =  —C(A  - 1)  .  Good  quantitative  agreement  is  found  in  Fig. 
2  with  respect  to  experimental  data  (given  in  Table  3).  The  relationship  (4),  with  v+  =  18.2  m/s 

and  c+  =  1,600  ml s provides  results  (dashed  line)  which  are  very  close  to  those  predicted  by 

using  the  expression  (39)  of  the  Hugoniot  extrapolated  to  large  stresses  (up  to  10  GPa).  This 
comparison  provides  a  first  validation  of  the  foregoing  extrapolation. 

5)  Structure  of  plastic  shocks  in  laminates 

For  stress  amplitude  cr  larger  than  1  GPa,  the  relationship  (41)  is  not  satisfied,  and  states 
(0)  and  (-)  are  related  by  a  single  shock  as  shown  in  Fig.  6,  composed  of  an  elastic  jump  from 
state  (0)  to  (*)  and  plastic  shock  from  state  (*)  to  (-).  The  value  of  the  stretch  A*  at  state  (*) 

corresponds  to  the  intercept  of  the  Rayleigh  line  (42)  with  the  elastic  path  A  —>  a  =  -T(A,1) 


emanating  from  state  (0). 


a\x  1)=-ru‘,i)  (44) 

For  the  amplitude  <j  =  2  MPa  ,  the  following  values  are  found,  £  =  0.981 
a*  =  0.28  MPa  . 

Consider  now  any  point  on  the  Rayleigh  line  between  states  (*)  and  (-)  i.e.,  for  X~  <  X  <  £  . 

For  each  point  there  is  a  one-to-one  correspondence  between  X  and  the  plastic  stretch  Xp  in 
the  interval  Xp  <XP<\ .  This  correspondence  is  given  by, 

Tx{X,Xp)  =  p0Cz(X- 1).  (45) 

which  can  be  resolved  into, 

A  =gN(Ap).  (46) 

Similarly,  on  the  Hugoniot,  there  is  a  one-to-one  correspondence  between  X  and  Xp  denoted 
as, 

A  =  gH(Xp)  (47) 

For  a  given  value  of  the  plastic  stretch  Xp ,  consider  the  intersections  ( I R )  and  ( IH  )  of  the 
elastic  path  X  — >  a  =  -T(X,XP)  with  the  Rayleigh  line  and  the  Hugoniot,  respectively,  see  Fig. 
7.  Note  that  the  same  value  of  the  plastic  stretch  Xp  is  associated  with  I H  and  I R ,  and  that 
&  =gH(AP) for  !h  and  X  =gR(Xp )  for  IR. 

The  compressive  stress  <jIh  = -T(gH(Xp),Xp)  associated  with  IH  is  the  yield  stress 

obtained  at  the  rear  of  the  plastic  front  in  a  plate  impact  experiment  with  stress  amplitude  <jIh  . 

Therefore,  this  yield  stress  can  be  associated  to  a  state  of  homogeneous  plastic  deformation 
defined  by  the  plateau  behind  the  plastic  shock.  However,  within  the  shock  structure  the  plastic 
deformation  is  highly  non-uniform.  For  these  high  plastic  deformation  gradients,  conventional 
plasticity  laws  are  not  able  to  restitute  properly  the  overall  material  response,  in  particular  when 
the  shock-width  is  comparable  to  the  characteristic  lengths  of  the  composite  internal  structure.  To 


account  for  the  interaction  between  the  internal  structure  and  the  propagating  wave,  it  is 
appealing  to  consider  a  nonlocal  plasticity  theory. 

Nonlocal  theories  of  plasticity  have  been  developed  in  the  last  two  decades  to  fonnulate 
macroscopic  phenomenological  laws  taking  into  account  the  effects  of  the  material  internal 
structure,  see  for  example,  Aifantis  (2003),  Fleck  and  Hutchinson  (1993,  2001),  Gao  et  al.  (1999), 
Huang  et  al.  (2000).  Our  aim  here  is  to  model  geometrical  dispersion  in  laminate  media  by  using 
nonlocal  plasticity  models.  Two  gradient  plasticity  fonnulations  and  their  predictions  are 
compared. 

Dissipation  of  mechanical  work  into  heat  by  viscosity  can  also  contribute  together  with 
dispersion  effects  to  the  build-up  of  the  shock  structure.  To  separate  these  effects,  viscous 
dissipation  will  not  be  considered  here. 

5,1)  First  order  gradient  plasticity  approach 

In  this  section,  a  first  order  gradient  plasticity  approach  is  considered.  Such  theories  have 
been  recently  developed  in  the  context  of  polycrystalline  materials,  the  internal  characteristic 

length  being  of  the  order  of  the  grain  size  (about  10 /urn  ),  Fleck  and  Hutchinson  (1993,2001),  Gao 

et  al.  (1999),  Huang  et  al  (2000).  For  the  laminated  media  considered  here,  the  internal  length  is 
expected  to  be  scaled  by  the  layer  widths. 

In  the  first  order  gradient  plasticity  theory  used  here,  the  yield  stress  <jy  is  a  function  of  the 

accumulated  plastic  strain  ep  and  of  the  first  gradient  of  this  strain.  It  is  assumed  that, 

<jy  =  aY0(sp)  +  K\psp  /  dx ^  (48) 

In  the  compression  tests  under  consideration  sp  can  be  identified  to  the  opposite  of  the 
longitudinal  plastic  strain. 

sp=-£p=- log(^)  (49) 

k  and  P  are  positive  material  parameters  which  are  related  to  the  internal  structure  of  the 

laminate  (geometry  of  the  layers,  material  characteristics  of  the  constituents,  including  the  mass 
densities  since  the  impedance  mismatch  between  layers  plays  a  role  in  wave  reflection  and 
refraction). 

The  yield  stress  crY  corresponding  to  the  state  IR  differs  from  the  yield  stress 


< 7y0  associated  with  the  state  IH  on  the  Hugoniot,  the  latter  being  related  to  a  state  of 
homogeneous  deformation.  In  the  present  modeling,  the  gap  between  these  yield  stresses  is 


attributed  to  the  high  level  of  heterogeneity  of  plastic  deformation  in  the  shock  front.  The  yield 
stress  is  assumed  to  be  affected  by  the  plastic  strain  gradient  as  in  (48), 

<jJr  =  <jIh  +  ic\d£p  / dx p  (50) 

with, 

<?Ih  =  -T(gH(Ap),Ap),  °iR  =-T(gR(*PUP).  (51) 

Defining  the  function, 

f(sr)  =  (-(-T(gR(r)X,)  +  T(gH(r\X")))'lf  (52) 

K 

where  Ap=exp(ep)=exp(-ep) ,  (50  )  can  be  rewritten  as, 

dsp 

=  ~f(ep ).  (53) 

ox 

The  sign  (-)  comes  from  the  fact  that  sp  is  a  decreasing  function  of  x  for  a  wave  moving  to  the 
right.  Note  that  the  plastic  deformation  varies  within  the  shock  from  sp*= 0  (state  *)  to 
sp ~  =  -  \og(Ap~)  (state  -)  and  that  f(ep )  vanishes  for  state  (-). 

Upon  integration  of  (53)  between  state  (*)  and  state  (-),  the  shock  profile  can  be  obtained  for  the 
particle  velocity,  the  total  strain  s  =  \og(A)  and  the  cumulated  plastic  strain  sp  =  -log(//' ) . 

All  numerical  calculations  in  the  following  are  made  for  a  linear  dependence  of  the  yield 
stress  upon  the  plastic  strain  gradient,  i.e.,  for  /?  =  1  ,  and  a  good  correspondence  with 

experimental  results  is  found  for  that  value  of [i .  However,  in  general  ft  can  have  any  positive 

value  and  the  effect  obtained  by  varying  the  value  of  ft  will  be  considered  later  when 
discussing  second  gradient  plasticity  approaches. 

Figure  8  shows  the  profile  of  sp  within  the  shock  and  the  variation  of  j[£p)  .  Calculation 

were  made  for  the  stress  amplitude  a  =  2  GPa  and  k  =  1 .7  m  MPa  ,  (5  =  1  .  A  smooth 

transition  is  observed  when  E,  — >  -qo  .  On  the  contrary  the  value  sp  =  0  corresponding  to  state 

(*)  is  reached  at  a  finite  value  of  E,  and  for  a  non-zero  value  of  the  slope  dspldE,. 

Note  that  according  to  the  present  modeling,  there  is  no  elastic  precursor.  This  is  in 


agreement  with  experimental  observations  of  Zhuang  et  al  (2003). 

An  important  role  in  structuring  the  shock  is  played  by  the  parameter  k  which  accounts  for 
the  internal  structure  of  the  laminate.  From  (48)  it  follows  that  k  has  the  dimension  of 

(length)^  x  (stress)  .  Considering  that  k  is  a  function  of  the  layer's  geometry  (specified  by  L, 

and  L=Ll+L2)  and  its  material  properties,  and  using  dimensional  analysis  (Vaschy-Buckingham 
theorem),  it  is  easily  found  that, 


k{Lp  orefyx  =  F(Ll  /  L,  non  -  dimensional  material  parameters)  (54) 

F  is  a  non-dimensional  quantity  and  alvf  is  any  characteristic  stress,  for  instance  the  elastic  shear 
modulus  of  phase  1.  Non-dimensional  material  parameters  are  of  the  form:  A1  /  <jref ,  px  /  <jref, 

A  /  crref ,  ju2  /  oref ,  (j n  I  a ref  ,  /;,  ,  crY2  /  crref  ,  n2 ,  pjp2,  where  Ai ,  pi  ,  uYi ,  /?,  are  for 

phase  (i)  respectively:  the  two  Lame  constants,  the  initial  yield  stress  and  the  hardening  exponent 
(  if  a  power  law  of  Hollomon  or  Ramberg-Osgood  type  is  assumed  to  describe  strain  hardening 

of  a  given  constituent).  Note  that  mass  densities  px  and  p2  of  constituents  are  involved  through 

the  effect  of  the  impedance  mismatch  in  reflection  and  refraction  of  waves  at  interfaces  beteen 

different  phases.  Mass  densities  are  combined  in  the  ratio  pj p2  ■  Therefore,  for  a  given  volume 

fraction  f=Ll/L  of  phase  1,  k  appears  to  be  proportional  to  Lp  based  on  (54). 


For  the  stress  amplitude  cr  =2  GPa  and  for  k  =  1 .7  m  Mpa  ,/?  =  !,  the  shock  width  can 


be  evaluated  approximately  as  0.5  mm  from  Fig.  8.  A  direct  estimate  of  the  shock  width  is 
obtained  from  the  relationship, 


£p- 

dsp  jdt; 

max 

(55) 


where  \sp\  is  the  maximum  plastic  strain  rate  within  the  shock  layer.  From  (53), 

I  I  max 

=  max  (f(£p )  .  Note  also  that  \sp\  =(j\dsp IdE I  .  For  <j  =2  GPa  the  shock 

I  I  max  0<£p<£p~  Imax  I  ^  I  max 

width  in  PC37/SS19  is  estimated  using  (55)  as  vv=0.49  mm  in  agreement  with  the  shock  profile 
shown  in  Fig.  8.  Note  that  w  is  comparable  to  the  cell  size  L= 0.56  mm. 

For  a  given  volume  fraction  of  constituents,  the  following  relationships  between  the  internal 
structure  of  the  composite  material  and  the  shock  structure  can  be  demonstrated: 


K  OC  Lp , 


(56) 


sp  oc-)- ,  wccL 

max  Jj 


The  proportionality  of  k  with  respect  to  if  results  from  (54).  By  using  (53)  one  obtains 


ids  i 


e  p=C^=-Cj{ s ,,)  ,  and  from  the  dependence  k  cc  if  ,  it  follows  with  (52)  that 


sp  oc  k  oc  — .  Finally,  with  (55),  the  width  of  the  shock  layer  appears  to  be  proportional  to 

max 


L,  irrespective  of  the  value  of  the  exponent  ft . 


In  Fig.  9  the  stress  amplitude  is  displayed  in  terms  of  sp  in  a  log-log  diagram,  for 


/?  =  1 .  Results  of  the  model,  given  by  the  solid  line  for  PC37/SS19  and  by  the  dashed  line  for 

PC74/SS37,  reproduce  quite  well  the  experimental  data.  Theoretical  predictions  are  obtained  by 
using  material  characteristics  given  in  Tables  2  and  4. 

PC37/SS19  and  PC74/SS37  laminates  have  same  proportion  of  constituents,  i.e.,  same 
volume  fractions.  The  only  difference  between  those  laminates  is  the  layer's  total  thickness 
which  is  double  for  PC74/SS37  in  comparison  to  PC37/SS19.  The  cell  size  has  the  value 

L  =  0.56  mm  for  PC37/SS 19  and  L  =  1.12  mm  for  PC74/SS37.  For  a  fixed  value  of  the  volume 


fraction/  and  for  /?  =  1 ,  k  is  proportional  to  L,  according  to  (54).  Therefore,  the  value 

k  =  3.4x1 06  m.Pa  for  PC74/SS37  is  taken  to  be  twice  the  value  k  =  1.7x1 06  m.Pa  for 

PC37/SS19.  Since  by  (56),  \sp\  oc l//rx  1/L ,  the  strain  rate  \ep\  for  a  given  stress  amplitude 

a-,  has  a  value  twice  larger  for  PC37/SS19  than  for  PC74/SS37.  This  is  in  agreement  with 
experimental  data  shown  in  Fig.  9. 

Results  in  Fig.  9  can  be  used  to  represent  both  laminates  through  the  following  power  law: 

\sp  I  (57) 

with  hL=  1.82  .  It  can  be  shown  that  the  slope  1  /hL  of  the  lines  in  Fig.  9  does  not  depend  upon 

k  but  is  a  function  of  (5  (this  will  be  further  discussed  when  considering  second  gradient 

plasticity  theories).  However,  for  a  given  value  of  /? ,  the  position  of  those  lines  does  depend  on 
k  .  Note  that  the  first  gradient  theory  with  a  linear  dependence  of  the  yield  stress  upon  the  plastic 


strain  gradient  (/?  =  1)  provides  a  good  prediction  of  \!h,  together  with  the  appropriate  effect  of 

the  cell  size  on  strain  rates.  The  value  of  k  for  the  PC37/SS19  laminate  was  calibrated  so  as  to 

get  the  best  fit  with  experimental  results  (re  =  1.7  m.MPa).  The  comparison  with  experimental 

results  for  PC74/SS37,  Fig.  9,  shows  a  quite  good  agreement  and  furnish  an  additional  check  of 
the  pertinence  of  the  present  modelling.  Indeed,  it  is  worth  mentioning  that  experimental  results 
in  Fig.  9,  obtained  for  two  laminates  subjected  to  various  shock  intensities,  have  been  described 

for  (3  =  1 ,  with  a  single  calibration  of  the  parameter  k  for  one  of  the  laminates,  the  value  of 

k  for  the  second  laminate  being  deduced  from  the  theory  (  k  =  3.4  m.MPa  for  PC74/SS37which 
is  twice  larger  than  for  PC37/SS19). 

The  total  strain  rate  within  the  shock  can  be  estimated  as  the  ratio  of  the  total  strain  variation 
across  the  shock  s  by  the  rise  time,  w/C , 

£=d\£-^W 

Using  the  estimate  (54)  of  the  shock  width, 

~  £  |  | 

£=— k 

gP  I  Imax 

In  Fig.  10  the  dependence  of  the  stress  amplitude  cr  upon  \ep\  and  s  is  shown 

in  a  log-log  diagram.  A  small  shift  towards  larger  strain  rate  is  obtained  when  considering 

s  .  However  the  value  hL= 1.82  ,  which  characterizes  the  slope,  is  essentially  left  unchanged. 

The  small  difference  between  the  results  of  Fig.  10,  presented  either  in  terms  of  the  total 
strain  rate  or  in  terms  of  the  plastic  strain  rate,  justify  that  no  distinction  is  made  in  the 
discussion  between  the  representations  (1)  and  (57). 

Figure  1 1  shows  the  evolution  of  the  shock  width  w  in  terms  of  the  stress  amplitude 

a  for  PC37/SS19  and  PC74/SS37  laminates.  Since  the  volume  fraction  of  constituents  is 

identical  for  the  two  laminates,  w  is  proportional  to  L  (0.56  mm  for  PC37/SS19,  twice 
larger  for  PC74/SS37).  Thus  for  a  given  stress  amplitude,  w  is  twice  larger  for  PC74/SS37. 

The  use  of  a  nonlocal  constitutive  law  such  as  the  one  proposed  in  (50)  is  of  interest 
when  the  characteristic  length  of  stress  (or  strain)  heterogeneity  (which  can  be  identified 

here  with  the  width  of  the  shock  layer  w  )  is  of  the  order  of  the  characteristic  lengths  L]  and 


(58) 

(59) 


L2  defining  the  laminate’s  internal  structure.  Considering  for  instance  the  PC37/SS19 


laminate,  values  of  w ( from  0.984  mm  at  a  =\GPa  to  0.091  mm  at  cr-=\0GPa)  are 


comparable  to  Lx  =  0.37  mm  and  L2  =0.19  mm . 

The  effect  of  second  order  elastic  constants  is  on  the  results  is  shown  in  Fig.  12(a)  where 
the  stress  amplitude  cr  is  reported  in  terms  of  the  maximum  of  the  plastic  strain  rate 

\£p\  (log-log  diagram).  In  the  elastic  law  (17),  second  order  elastic  constants  are  scaled  by 

the  parameter  b.  Previous  results  obtained  for  the  value  6=66. 7  (solid  line)  are  compared 
with  those  for  b=\  00  (dashed  line).  A  little  shift  and  a  small  change  in  the  slope  are  obtained 
between  the  two  lines  representing  results  for  b=66.7  and  b=100.  Similar  trends  were 
obtained  when  analyzing  shock  waves  in  polycrystalline  metals,  see  Fig.  15  in  Molinari  and 
Ravichandran  (2004). 

The  effect  of  /t+  (defining  the  state  (+))  is  shown  in  Fig.  12(b).  The  solid  line 
corresponds  to  the  previous  results  of  Fig.  9  (a  versus  \£p\  with  A'  =0.99 ).  Two 
different  curves  are  shown  for  /t+=0.97 .  The  solid  line  with  dots  represents  cr-  in  terms  of 

\sp\  The  dashed  line  shows  the  variation  of  the  stress  jump  Acr  =  a  -  cr+  with  respect 

I  Imax 

to  \sp\  .  This  result  indicates  that  the  effect  of  X  is  almost  negligible  when  considering 
variations  of  Acr  in  terms  of  \£p\  .  Note  that  for  X=0. 99 ,  cr+ =0.1  \5GPa  is  negligible 

with  respect  to  a  for  the  range  of  stress  amplitudes  considered  here;  thus  cr  -Act  . 

However,  this  is  not  the  case  for  /t+=0.97  ,  for  which  cr+=  0.63  lGPa  .  From  this 

discussion,  it  can  be  concluded  that  results  presented  in  Fig.  9  would  only  be  slightly 
affected  by  considering  values  of  X  smaller  than  0.99. 


5.2)  Second  order  gradient  plasticity  approach 

Second  order  gradient  plasticity  approaches  are  frequently  used  in  the  literature  to  analyze 
strain  localization  phenomena,  for  example,  Aifantis  (2003).  In  this  section  the  pertinence  of 
such  constitutive  modeling  with  respect  to  shock  wave  analysis  in  laminate  media  is  evaluated. 

The  yield  stress  ay  which  was  given  by  (48)  for  the  first  order  gradient  approach  is 
presently  assumed  to  be  dependent  upon  the  second  gradient  of  the  accumulated  plastic  strain, 


o>  =(TY0(£  )  +  /cl 


2rP 


ds 


dx 2 


(60) 


where  kx  and  /?,  are  positive  material  parameters  accounting  for  the  laminate  internal  structure. 
The  relationships  (50)  and  (53)  are  replaced  by. 


<j,  =  cr,  +  k, 

2R  1H  1 


2  7 ~P 


ds 


dx2 


(61) 


and 


d2s 1 
dx2 


-=f  (sp) 


(62) 


with 


l  (?')=(— {.-ngMr\w)+ngH(w\v)))'lf' . 

K, 


(63) 


dsp 

First  quadrature  of  (62)  can  be  obtained  by  multiplying  both  sides  by  - .  Upon  integration 

dx 


from  state  (*)  to  the  present  state,  and  assuming  a  smooth  transition  at  state  (*),  i.e., 


ds1 


dx 


=0, 


dsp 

and  a  negative  slope  - <  0  elsewhere  in  the  shock  layer,  one  obtains, 

OX 


d£P 

dx 


-V2£(o 


(64) 


with 


g(£P)=\h  (ddg 


(65) 


For  the  stress  amplitude  cr  =  2  GPa  ,  Fig.  13  shows  how  shock  profiles  compare  for  the 
first  gradient  and  second  gradient  plasticity  theories,  defined  by  (48)  and  (60), respectively  with 
parameters  ( k  =  1.7  m.MPa  ,/?  =  !)  and  (  kx  =  750  m2. Pa  ,  /?,  =  1 ).  For  the  second  gradient  model, 


a  smooth  transition  to  sp  =  0  is  observed  when  E,  increases.  On  the  contrary  the  value 

sp =0. 118  corresponding  to  state  (-)  is  reached  at  a  finite  value  of  E,  beyond  which  an 
homogeneous  state  is  obtained.  The  reverse  situation  is  observed  for  the  first  gradient  theory;  the 
state  (-)  being  reached  asymptotically  when  E,— >  -qo  ,  while  the  state  (*)  corresponds  to  a  finite 

value  of  E,  with  a  slope  different  from  zero. 

In  Fig.  14  the  profiles  of  plastic  strain  ep  and  the  magnitude  of  total  strain  |^|=|log(A)|  are 
shown  for  the  second  gradient  approach  ( kx  =  750  nr. Pa  ,  f\  =1).  The  stress  amplitude  is 

a  =  2  GPa  .  At  the  shock  front,  there  is  an  elastic  jump  from  state  (0)  to  state  (*),  see  Fig.  6. 

Thus,  the  total  strain  is  discontinuous  at  the  shock  front,  while  the  plastic  strain  is  continuous  and 
equal  to  zero. 

The  parameter  kx  in  (60),  accounting  for  the  effect  of  the  laminate’s  internal  structure  on 
the  constitutive  response,  follows  the  following  relationship  obtained  from  dimensional  analysis, 
kx  =  F(L]  /  L,  non  -  dimensional  material  parameters)  (66) 

It  appears  that,  for  a  given  volume  fraction  of  constituents,  kx  is  proportional  to  Z.2/'  .  Note  that 

for  the  first  gradient  model,  k  was  found  to  be  proportional  to  Z/; . 

The  shock  width  wx  can  be  estimated  by  using  (55)  and  (64).  Thus,  it  is  found  that, 


and 


d£p 

dE, 

max 

\sp-\ 

d£p/d% 

max 

V2g(F~). 


(67) 


(68) 


For  the  PC37/SS 19  laminate  with  kx  =750  m2.Pa  and  J3X  =  1 ,  the  shock  width  w;=0. 385 


mm  is  obtained  from  (68)  for  a  =2  GPa  . 


For  a  given  volume  fraction  of  constituents,  it  follows  immediately  from  the  definition  (65) 
of  g  and  from  the  dependence  k  oc  L1,:  that, 

w{ccL,  \£p\  cc-j-  (69) 

I  Imax 

The  scaling  of  the  shock  width  and  of  the  magnitude  of  plastic  strain  rate  by  the  laminate 
cell  size  L  ,  is  similar  to  the  scaling  obtained  in  (56)  with  the  first  gradient  theory. 

The  major  difference  between  the  two  gradient  plasticity  approaches  appears  when  looking 

at  the  dependence  of  the  stress  amplitude  cj  upon  \sp\  shown  Fig.  15.  Compared  to 

experimental  results  of  PC37/SS19,  it  appears  that  the  slope  1  /hL  is  not  adequate  for  the  second 

gradient  model  with  /?,  =  1  (linear  dependence  of  the  yield  stress  upon  the  second  gradient  of 
the  plastic  deformation). 

Figure  16  illustrates  the  effects  of  parameters  /?,  and  kx  on  the  relationship  between  stress 

amplitude  and  magnitude  of  the  plastic  strain  rate.  For  a  given  value  of  (/?,  ,kx)  this  relationship 
is  fairly  well  represented  in  the  log-log  diagram  of  Fig.  16  by  a  straight  line  denoted  by 
L(/?,  ,kx  )  with  slope  1 1 hL  .  Calculations  show  that  this  slope  depends  on  the  value  of  /?,  but  is 

nearly  independent  of  kx  ,  see  lines  A,  B  and  C,  D  in  Fig.  16.  For  /?,  fixed,  L(/?,  ,kx)  is  shifted 

horizontally  by  changing  the  value  of  kx  .  Thus,  there  is  a  one-to-one  correspondence  between 

couples  (/?,  ,kx  )  (/?,  >  0,/fj  >0)  and  straight  lines  L(/?,  ,kx  )of  positive  slope  in  the  log-log 

diagram  of  Fig.  16.  The  material  parameters  ( ,kx)  defining  nonlocal  plasticity  can  be 

uniquely  defined  from  the  best  fit  of  L(/?j  ,  kx  )  with  experimental  data. 

From  (66)  and  (69)  it  follows  that,  for  given  value  of/?,  and  of  the  stress  amplitude  <j  , 

\£p\  cc  k \'2Pi  .  Thus,  for  two  different  values kxa  and  kxb  the  magnitudes  of  plastic  strain  rate 

corresponding  to  a  given  cr  are  related  by  the  following  relationship, 

-1/2  a 

(70) 


Note  that  this  quantity  is  not  dependent  on  the  stress  amplitude  cr  . 

Therefore,  considering  the  lines  A  and  B  in  Fig.  16,  respectively,  associated  with 
(/?,  =1  ,kxa  =  103)  and  (/ 31  =1  ,kxb  =412)  ,  one  has  \£PA\  =  2£f  for  any  value  of  a  . 

Similarly,  between  lines  C  (y 3X  =0.5  ,kxc  =  7.5x105A)  and  D  (/?j  =0.5  ,kw  =  15x10^A)  one 

has  \eZ\  =  2|fd 

I  ^  I  max  I  u  I  max 

For  a  given  characteristic  length  L*  and  stress  Z  ,  consider  the  family  of  lines  defined  as, 

F  (L* ,  E)  =  {  L(/?j ,  kx  )  such  that  E  =  k \  l8~2p'  }  (71) 

It  can  be  motivated  from  scaling  arguments  and  verified  numerically  that  all  the  lines  of 
F  (L*,Z)  have  an  unique  intersection  point.  Let  us  denote  by  F,  and  F;  the  families  with 

intersection  points  respectively  I  and  J.  Examples  are  shown  in  Fig.  16  of  families  F,  (  F,  )  with 
lines  A  and  C  converging  at  I  (lines  D  and  B  converging  at  J).  In  the  example  shown  in  Fig.  16 
the  characteristic  length  associated  to  F:  has  the  value  L*  =0.137  mm  and  the  value  of  L*A  is 
double  that  of  the  above  value  (0.274  mm). 

Consider  in  general  the  family  F,  obtained  by  horizontal  translation  of  F,  with 
magnitude  defined  by  8  =  ef  I  sp ,  as  illustrated  in  Fig.  16  for  8  =  2.  Defining  by  (L^Z,) 

(  (Ljj,Ey))  the  characteristic  length  and  stress  associated  to  F,  (  F,  ),  it  is  proven  in  Appendix 
A  thatZ*  / L]  =  8  and  Z7  =  Z7  . 

Similar  considerations  could  be  developed  for  the  first  gradient  plasticity  model  (48).  A 
one-to-one  correspondence  can  also  be  found  between  the  parameters  ( /?  ,  k  )  defining 
non-local  plasticity  and  straight  lines  with  positive  slope  in  the  diagrams  of  Fig.  15  or  Fig.  16. 
Figure  17  shows  the  variation  of  the  stress  amplitude  cr  in  terms  of  the  maximum  of  the 

plastic  strain  rate  for  the  two  laminates  PC37/SS19  and  PC74/SS37.  Experimental  results  are 
well  described  by  the  second  gradient  plasticity  model  (60)  with  /?,  =  0.5  and  kx  =7.5x10:>  N 


for  PC37/SS19,  and  kx  =  15x105N  for  PC74/SS37.  As  for  the  first  gradient  model,  atx  is 


scaled  by  the  cell  size  (the  value  of  kx  is  doubled  for  PC74/SS37). 

Results  obtained  for  the  PC37/SS19  laminate  are  summarized  in  Table  5  for  the  two  gradient 
plasticity  theories,  and  different  values  of  the  structure  parameters  ft  and  k  .  It  is  shown  in  Fig. 
18  (see  also  Table  5)  that  values  of  structure  parameters  can  be  chosen  so  as  to  provide  almost 
identical  results  for  the  first  gradient  ( /?  =  0.88  ,  k  =  3.55  MPa)  and  the  second  gradient 

approaches  ( /3X  =  0.5 ,  /y  =  7.5x10'’  N).  Therefore,  the  relationship  between  stress  amplitude 

and  maximum  plastic  strain  rate  is  enough  to  characterize  the  structure  parameters  for  a  given 
gradient  plasticity  approach  but  is  not  sufficient  to  discriminate  between  the  two  plasticity 
theories  considered  here.  The  shape  of  the  shock  front  provides  additional  information  allowing 
favoring  one  of  the  theories.  The  smooth  transition  towards  state  (-)  shown  in  Fig.  13  for  the  first 
gradient  theory  seems  to  be  in  better  agreement  with  experimental  shock  profiles,  Zhuang  (2002, 
2003). 

6)  Summary  and  conclusions 

First  and  second  gradient  plasticity  theories  were  used  to  describe  steady  shock  waves 
propagating  through  laminates  in  the  direction  perpendicular  to  layers.  The  aim  of  this 
investigation  was  to  analyze  and  model  experimental  results  obtained  for  two  laminates  made  of 
the  same  constituents.  The  only  difference  being  that  the  layer’s  widths  in  the  laminate 
PC74/SS37  are  doubled  when  compared  to  those  of  PC37/SS19.  The  maximum  strain  rate  within 
the  shock  layer  was  measured  to  be  twice  smaller  for  PC74/SS37.  Experimental  results  could  be 

represented  by  the  power  law  (57),  \ep\  =  B[a  r  where  I sp\  is  the  maximum  plastic 

strain  rate  and  cr  is  the  stress  amplitude.  The  exponent  h,  has  the  same  value  for  both 

laminates  (close  to  1.8)  ,  however  B  is  twice  larger  for  PC37/SS19  in  comparion  to  the  value  for 
PC74/SS37. 

These  experimental  results  could  be  well  described  by  using  the  first  gradient  plasticity 

model  (48).  In  this  approach,  structure  parameters  k  and  ft  are  introduced  to  account  for  the 

internal  structure  of  the  laminate  (layer's  thickness  and  material  properties  of  constituents). 

Since  the  volume  fractions  of  phases  are  identical  for  both  laminates,  the  difference  in  the 
material  response  can  be  related  to  the  laminate  cell  size,  L  .  The  present  modeling  shows 

thatl^l  cc-j-  .  The  experimental  finding  that  the  strain  rate  within  the  shock  layer  is  twice  larger 
in  PC37/SS19  than  in  PC74/SS37  is  then  rationalized. 


Similar  results  were  found  by  using  the  second  order  gradient  plasticity  approach  (60). 
Experimental  results  are  well  described  by  taking  a  non-linear  dependence  of  the  flow  stress  with 

respect  to  the  second  gradient  of  plastic  deformation,  with  exponent  (3X  =0.5  in  (60). 

Non-local  effects  are  described  in  these  gradient  plasticity  theories  by  two  material 

parameters,  (3  and  k  .  It  was  shown,  that  for  a  given  plasticity  theory,  these  parameters  ( (3 .  k  ) 

could  be  uniquely  detennined  from  experimental  results.  However,  both  plasticity  theories 
provide  results  of  comparable  quality  when  compared  to  experimental  data  relating  the  stress 
amplitude  to  the  maximum  of  the  plastic  strain  rate.  Nevertheless,  the  first  gradient  plasticity 
approach  can  be  favored  when  comparing  the  shape  of  the  shock  front  to  experimental  data. 

Viscous  effects,  which  have  been  neglected  in  the  present  modeling,  have  sometimes  a  major 
role  in  structuring  shock  fronts  as  for  example  in  polycrystalline  metals.  In  the  latter  case, 
experimental  results  show  a  power  law  dependence  between  the  stress  jump  across  the  plastic 
shock  and  the  maximum  strain  rate,  with  exponent  h= 4 ,  Swegle  and  Grady  (1985),  a  value 

much  larger  than  those  of  hL  for  laminates  (of  the  order  of  2). 

The  present  theory  can  be  augmented  to  account  for  viscous  effects.  However,  the 
experiments  performed  on  laminates  indicate  that  the  scaling  of  the  shock  layer  is  mostly  related 
to  the  laminate  internal  structure  (in  particular  to  the  cell  size).  Viscous  effects,  which,  for  metals 
and  polymers,  can  be  related  to  basic  mechanisms  of  visco-plastic  deformation,  are  not  able  to 
account  for  such  scaling.  Dispersion  effects  due  to  heterogeneities  of  the  laminate  structure  are 
thought  to  be  all  the  more  important,  as  the  contrast  between  the  mechanical  properties  of 
constituents  is  larger.  In  the  case  of  a  large  contrast,  wave  dispersion  appears  to  be  the  dominant 
factor  structuring  the  shock  front  and  viscosity  is  believed  to  play  a  minor  role. 


Table  1 :  Material  characteristics  of  polycarbonate  and  stainless  steel 


Poisson  ratio 

Elastic  shear  modulus 

Mass  density 

PC  (Polycarbonate) 

0.37 

0.94  Gpa 

11 90  kg/m3 

SS  (Stainless  steel) 

0.29 

77  Gpa 

7890  kg/m3 

Table  2:  Overall  material  characteristics  of  the  laminates  PC37/SS19  and  PC74/SS37 

Gin 

G 122 

b 

Po 

5  c+ 

PC37/SS19  and  PC74/SS37  6.84  Gpa 

3.60  Gpa 

66.7 

3468  kg/m3 

2.1  1600m/s 

Table  3:  Experimental  data  for  the  shock  velocity  C  and  the  particle  velocity  v- 


PC37/SS19 

C  (km/s)  1.96 

1.89 

2.39 

2.32 

2.33 

2.53 

3.12 

n  /  \  0.254 
v-  (km/s) 

0.254 

0.454 

0.454 

0.485 

0.485 

0.699 

PC74/SS37 

C  (km/s)  1.94 

2.39 

n  /  \  0.257 
v-  (km/s) 

0.527 

Table  4:  Parameters  related  to  the  laminate  internal  structure 


A 

^2 

K 

PC37/SS19 

0.37  mm 

0.19  mm 

1.7 

rn.MPa 

PC74/SS37 

0.74  mm 

0.37  mm 

3.4 

rn.MPa 

Table  5:  Theoretical  results  for  the  PC37/SS19  laminate  using  the  values  of  material  parameters  given  in 


Table  2. 


First  gradient 

k=\  .7  m  MPa 

Second  gradient 

Ay  =750  m2  MPa 

a(GPa ) 

v-(m/s ) 

C(m/s) 

\ep\  (s'xlO6) 

1  Imax 

w  (mm) 

AH  (s'xlO6) 

1  Imax 

wl  (mm) 

8 

748 

3084 

6.15 

0.117 

3.93 

0.183 

6 

615 

2811 

3.69 

0.160 

2.75 

0.215 

4 

461 

2497 

1.79 

0.246 

1.64 

0.268 

2 

273 

2112 

0.51 

0.490 

0.648 

0.385 

Impact  direction 


L 


2 


Fig.  1 :  Schematic  of  the  geometry  of  a  periodically  layered  composite  with  layer's  thicknesses  L ,  and  L2 


Fig. 2:  Shock  velocity  C  versus  particle  velocity  v~  at  the  rear  of  the  plastic  shock.  Experimental  data  are 
represented  by  stars  and  squares  respectively  for  PC37/SS19  and  PC74//SS37  laminates.  The  dashed  line  is 
given  by  equ.  (4)  with  slope  <S=2.1,  v+=18.2 m/s  and  c+=l600m/s  .  The  solid  line  with  dots 


corresponds  to  calculations  made  in  section  4.2  by  using  the  Hugoniot  defined  by  (39) 


Fig. 3:  Decomposition  of  the  deformation  gradient  F  into  plastic  part  F_  P  and  elastic  part  F_' 


COMPRESSIVE  STRESS  a  (Gpa) 


Fig. 4:  Plastic  shock  for  the  stress  amplitude  a  =0.5GPa  .  The  material  is  first  compressed  from  the  initial 


state  to  the  state  (+)  by  an  elastic  shock  (elastic  precursor),  and  is  further  compressed  by  a  plastic  shock  from 
state  (+)  to  state  (-).  These  two  distinct  shocks  exist  if  the  condition  (41)  on  shock  velocities  is  satisfied. 


a  (Gpa) 


Fig.5:  Hugoniot  curve  (solid  line)  and  elastic  loading  path  A—>cr=— Tn(A,Y)  (dashed  line)  for  the  laminate 

PC37/SS 19.  Results  are  identical  for  PC74/SS37.  A  zoom  is  shown  for  0.97<A<1  in  the  upper  part  of  the 
figure. 


COMPRESSIVE  STRESS  a  (Gpa) 


Fig. 6:  Elasto-plastic  shock  generated  for  a  large  enough  stress  amplitude  such  that  condition  (41)  is  not 
satisfied. 


Fig. 7:  The  elastic  loading  path  corresponding  to  the  fixed  value  of  the  plastic  stretch  Ap  =0.96  intercepts  the 
Flugoniot  at  IH  and  the  Rayleigh  line  at  IR  .  The  shock  stress  amplitude  is  <j  =2GPa  . 


Fig.8:  Structure  of  the  plastic  shock  for  a  =  2GPa  : 

(a)  distribution  of  the  cumulated  plastic  deformation  sp  in  terms  of  the  position  q 

(b)  variation  of  the  function  f[sp)  defined  by  (52) 


The  first  order  gradient  plasticity  theory  is  used  here  with  K  =  1.7  m  Mpa  ,  /?  =  I 


Fig. 9:  Stress  amplitude  cr  in  terms  of  the  maximum  of  the  plastic  strain  rate  within  the  shock  layer. 

Experimental  results  for  the  laminates  PC37/SS19  and  PC74/SS37  are  respectively  represented  by  stars  and 

squares  .  Results  of  the  first  gradient  modelling  are  obtained  for  (3  =  1 ;  K  =  1.7  m  Mpa  is  used  for 

PC37/SS19  and  tc=3.4m.MPa  for  PC74/SS37.  Note  that  doubling  the  value  of  K  for  PC74/SS37  is  the 
direct  consequence  of  the  layer's  thicknesses  of  PC74/SS37  being  twice  those  of  PC37/SS19. 


FiglO:  Stress  amplitude  <j  in  terms  respectively  of : 

-  the  maximum  of  the  plastic  strain  rate  \£p\  within  the  shock  layer 

I  Imax 

-  an  estimate  of  the  maximum  of  the  strain  rate  within  the  shock  layer  £  . 

Stars  represent  experimental  data  for  the  laminates  PC37/SS19.  Results  of  the  modelling  are  obtained  for  the 
first  gradient  theory  with  K=\.lm.MPa  and  /?  =  !. 


Fig.  1 1 :  Shock  width  w  estimated  with  (55).  Note  that  the  shock  width  for  PC74/SS37  is  twice  that  for 
PC37/SS19.  Dots  correspond  to  the  following  stress  amplitudes  (Gpa):  1,  1.5,  2,  3,  4,  6,  8,  9,  10. 


STRESS  o'  (GPa) 


Fig.  12:  Effects  of  second  order  elastic  constants  (a)  and  of  A+  (b) 

(a)  Stress  amplitude  cr  in  terms  of  the  plastic  strain  rate  \£p\  for  two  values  of  the  parameter  b 

I  I  max 

which  scales  second  order  elastic  constants.  Results  with  b=100  are  compared  with  the  previous 
results  of  Fig.9  obtained  with  b=66.7.  Stars  are  experimental  data  for  the  PC37/SS19  laminate. 

(b)  Effect  of  /J  (defining  the  state  (+)).  The  solid  line  corresponds  to  the  previous  results  of  Fig.9 
(cr-  versus  \sp\  with  A'  =0.99  ).  Two  different  curves  are  shown  for  A+  =0.97 .  The  solid 

I  I  max 

line  with  dots  represents  cr  in  terms  of  \£p\  The  dashed  line  shows  the  variation  of  the 

I  I  max 

stress  jump  Acr=cr  —  cr+  with  respect  to  \£p\ 

I  Imax 


i;(mm ) 


Fig.  13:  Distribution  within  the  shock  structure  of  the  cumulated  plastic  deformation  £p  in  terms  of  the 
position.  The  first  gradient  plasticity  theory  defined  by  (48)  with  ( (3  =  1 ,  K= 1 ,1m. MPa  )  is  compared  to  the 
second  gradient  theory  (60)  with  ( /?,  =  \,KX  =  150N ).  The  stress  amplitude  is  cr  =  2  GPa  . 


Fig.  14:  Shock  wave  profiles  calculated  for  <j  =2 GPa  with  the  second  gradient  approach 
(/? j  =  1 ,  ATj  =  750N ).  The  dashed  and  the  solid  lines  correspond  respectively  to  distributions  of  the  plastic 
strain  sp  and  of  the  absolute  value  of  the  total  strain  |fi]=|log(/l)| . 


STRESS  o'  (GPa) 


Fig.  15:  Compared  to  the  first  gradient  model,  the  slope  of  the  (dashed)  line  associated  to  the  second  gradient 
model  (with  J3l  =  1 )  is  not  adequate  to  account  for  experimental  results  of  PC37/SS 19  (stars) 


MAXIMUM  PLASTIC  STRAIN  RATE  (x  1 06  s'1 ) 


Fig.  16:  Second  gradient  approach:  effects  of  parameters  and  Kx  on  the  relationship  between  stress 
amplitude  and  magnitude  of  plastic  strain  rate.  The  slope  depends  on  the  value  of  /3 \  but  is  not  dependent 


upon  ATj  .  For  a  given  value  of  ft  x  ,  the  lines  are  shifted  horizontally  by  changing  the  value  of  K  x  . 


Fig.  17  :  Stress  amplitude  <J  in  terms  of  the  maximum  of  the  plastic  strain  rate  within  the  shock  layer. 
Experimental  results  for  the  laminates  PC37/SS19  and  PC74/SS37  are  respectively  represented  by  stars  and 
squares.  Results  of  the  second  gradient  modelling  (60)  are  obtained  for  f5x  =0.5  and  Kx  =  7.5x105iV 

(PC37/SS19),  Kx  -  15x10s  N  (PC74/SS37).  As  for  the  first  gradient  model,  Kx  is  scaled  by  the  cell  size 
which  is  multiplied  by  two  for  PC74/SS37. 


STRESS  a  (GPa) 


Fig.  1 8  :  First  gradient  approach  (  /?  =  0.88  ,  K  =  3.55  MPa)  compared  to  second  gradient  approach 
( J3 j  =0.5,  Kx  =  7.5x1 0  N),  respectively  dashed  and  solid  lines.  Stars  represent  experimental  data  for 


PC37/SS19. 
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Appendix  A 

The  family  F  (L*j,Zy)  of  lines  intersecting  at  point  I,  is  translated  horizontally  as  shown  in 
Fig.  16  into  the  family  F  (Lj,Ey)  of  lines  intersecting  at  point  J.  Considering  in  F  (L*I,S/) 
the  lines  l_(/?1A,/ru)and  L (/?1C ,  at1c)  ,  denoted  for  simplicity  as  A  and  C,  one  obtains,  from  (71), 


=  kiaLj 


;-2 Aa  -  f- 2Ac 

—  Kx  CL,J 


(Al) 


The  lines  B  and  D  are  obtained  by  horizontal  translation  of  A  and  C  and  intersect  at  J, 
therefore  they  are  in  F  (L*, ,  X  y )  and  from  (71), 


2/  =  KIBLJ 


2  A, 


K\D^I 


'2  Ac 


where  the  following  relationship  has  been  used, 
P\  A  —  Pi B  1  PlC  —  P\D  ■ 


(A2) 

(A3) 


The  magnitude  of  the  translation  F  (L*j,Xy)— >  F  ( Lj,  X  y )  is  characterized  by  the  ratio 


a  =  £j  /  s j  .  Since  lines  A  and  C  intersect  at  /  and  B  and  D  intersect  at  J,  for  any  stress  amplitude 


cr  one  obtains, 


£P  £P 

°  A  °  r 

_ max  _  _ max  _  ^ 

kl  kl 

I  I  max  I  I  max 

From  (70),  it  follows  that, 


(A4) 


(A5) 


(A6) 


L]  =  aL) 

and  from  (Al)  and  (A2), 


(A7) 


(A8) 


